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Abstract: We construct a multi- field inflationary model consisting of multiple Kahler 
qj moduli derived from type IIB string compactification in the large volume limit. The model 

i consists of both heavy and light fields, with the former being frozen during the inflationary 

period and the latter acting as the inflaton(s). We study the evolution of all the fields 
during and after inflation until the preheating era when all the fields oscillate around their 
vacuum expectation values. Our numerical analysis shows that the curvature perturbations 
have an almost scale invariant power spectrum with n s ~ 0.96. 
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1. Introduction 

Ever since inflation was proposed in the early 1980s [1], there have been numerous infla- 
tionary models, including scalar field slow-roll inflation [2, 3]. (For a review, see [4, 5].) In 
recent years, there have been some promising developments in models derived or inspired 
by string theory, where scalar fields associated to the shape and size of the internal space, 
or to the positions of branes, serve as candidate inflaton fields. Based on the work of 
KKLT [6] and the large volume limit on moduli stabilization [7] [8] , an interesting slow-roll 
inflationary model has been proposed by Conlon and Quevedo [9], in which the inflaton 
is chosen to be one of the Kahler moduli, which sets the size of the compactified space. 
(For some more recent work along these lines, see [10] and [11].) Because the moduli only 
appear exponentially in the scalar potential, the potential is very flat along the directions 
of the light fields, which is ideal for obtaining successful inflation and a graceful exit to the 
normal Friedmann-Robertson- Walker (FRW) universe. 
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In this type of scenario, due to the presence of more than one modulus, we should 
expect to see a multi-field inflationary scenario in which all the fields play a role. In 
addition, due to the nature of cosmological models arising from string compactifications, 
the fields are not canonically normalized, i.e., the metric is in general neither diagonal nor 
field independent, and the scalar potential is highly non-trivial. This typical feature of 
string theory, leads to a highly coupled dynamical situation, in which, in principle, the 
motion of any one field impacts the evolution of the other fields. Therefore, we need to 
understand how each modulus and its perturbations evolve during and after inflation. In 
particular, we want to understand the difference between the "heavy" moduli and the 
"light" moduli. (For some earlier work, addressing the issue of the "very high energy 
physics" in inflation, see [12].) To accomplish this, it is necessary to explicitly solve the 
equations of motion for all the moduli and their perturbations. We study the resulting 
power spectra and the spectral indices of the perturbations, and compare the result with 
the usual single field scenario. 

In [9] the general idea of Kahler moduli inflation was outlined and a single modulus 
inflationary case was demonstrated by assuming that all other moduli are already in the 
final vacuum states. The same assumption was also made in [10] with two light fields. We 
will show that this assumption is redundant and may be dropped when all the moduli are 
taken into account. This is consistent with the numerical analysis in [11] which further 
pointed out that there exists a large region of parameter space within which the inflationary 
solutions fit the observations. Our work considers the general case of both multiple heavy 
and light fields. Due to the multi-field nature of our inflationary scenario, i.e., the existence 
of several light fields, we also compute the isocurvature perturbations. 

The outline of the paper is as follows. We first solve the perturbation equations of 
the scalar fields in Section 2. We also demonstrate the difference between the perturbation 
solutions of the heavy and light fields. In Section 3, we study the curvature and isocurvature 
perturbations. In Section 4, we introduce the moduli stabilization in type IIB string theory 
and the effective potential. In the numerical analysis, in Section 5, we investigate a few 
examples of the multi-field inflationary model and compute the spectra and tilts (spectral 
indices) which are used to make contact with observations. We conclude and summarize 
our results in Section 6. The detailed analysis of the Hankel functions, needed to study 
the effects of the heavy fields, can be found in Appendix A. 

2. Scalar Field Perturbations 

In this section, we study the scalar perturbations in a general multi-field model, following 
previous work, in particular that by Byrnes and Wands [13]. We explicitly show how the 
perturbations for the heavy fields are suppressed. 

2.1 The Background Equations of Motion 

Let us consider the scalar-tensor field theory of gravity, with the action of the standard 



form 
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where h ab is the metric on the space of fields, i^ a -space. From now on we will work in units 
where k 2 = 8irG N = M p ~ 2 = l 1 . 



The background spacetime metric g^ v is chosen to be the Friedmann- Robertson- Walker 
(FRW) metric. The line element is given by 



ds 2 



-dt 2 + a 2 (t)-fijdx' l dx j = a 2 {-q){-drf + jijdx l dx j ) 



where 



7y = 5ij[l + Akr 2 }- 2 . 
We will work in the flat universe in which k = 

( -1 \ 

a 2 {t) 

a 2 {t) 

\ a 2 {t) j 



(2.2) 
(2.3) 



(2.4) 



The background fields, g^ v and 4> a , satisfy the usual equations of motion obtained by 
varying the action S (2.1). This gives rise to the following field equations for the background 
fields (j) a , 



+ -^(d^g)\h ah g» v d v <p + V a = , 

dx^ y/—g 



(2.5) 



where -J^ is the covariant derivative. For flat space, g = —a 6 , and assuming that the fields 
are homogeneous, <j) a = 4> a (t), Eq. (2.5) becomes 



d 2 r , , # 6 # c ,, ab dv 



o, 



(2.6) 



where 7? is the Christoffel connection on the space of fields. By introducing the notation 
a = ifr> ^ a = d\^ a = < i )b ^' b4> a ? t ne field equations for <j) a can be further simplified 



4> a + 3H<p a + V' a = 0. 
Varying S with respect to the spacetime metric gives the Einstein equations 



3H 2 



H 



1 



+ V, 



(2.7) 

(2.8) 
(2.9) 



The background equations determine the evolution of the inflaton and how long infla- 
tion will last, the e-folding time, N. In most cases, the light fields and the heavy fields 3 are 



1 M P is the Planck mass, ~ 2.4 x 10 18 Gev. 

2 ^ denotes the covariant derivative on the coordinate space, while is the covariant derivative on the 
field space (labeled by b). 

3 The "light" fields have effective masses less than the Hubble parameter, m\ — Vn < V ~ H 2 , while 
the "heavy" fields are heavier than the Hubble parameter. 
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decoupled, i.e., the heavy fields remain constant most of the time when the light fields is 
rolling. Thus, the inflaton is always one, or a subset, of the light fields which controls the 
dynamics of inflation. However, in what follows, we will not make any assumptions about 
which fields are frozen and which are dynamic, but rather establish this fact as part of the 
calculation. 

2.2 Perturbations 

The most general scalar perturbations on the background spacetime metric is [13] [14] 

ds 2 = -(1 + 2A)dt 2 + 2aB- i dx i dt + a 2 [(l - 2ip) + 2E- ij ]dx i dx j (2.10) 
The associated equation of motion for the perturbed fields are [15] 

+ 3Hj t 5<p a + R a M H d 4> h + ^ a + W ah = ^J^^ ( 2 - n ) 

where the covariant derivatives and the Riemann curvature tensor are all evaluated on the 
field space. Here the fields, (f) a and their perturbations, 5<p a , are evaluated at a particular 
(comoving) wavenumber k, i.e., 

4> a = ^{t) = j d 3 xcj) a (x,t)e lkx (2.12) 

<{f(x,t) = J^XsJ d 3 kcf>t(t)e-* k * , (2.13) 

and similar for the perturbations. The perturbed Einstein equations are 

- 6H 2 A - 2k 2 HE = -A4> a 4> a + 4> a j^4>a (2.14) 
2HA = <p a 5<j) a . (2.15) 

To simplify the problem of solving the equations of motion for the perturbations, we 
work with the canonical field-space metric 4 in the spatially flat gauge. Eq. (2.11) then 
becomes [13] 

+ + pcP 1 + X>/ - V4j)W = • (2-16) 

To solve this equation, we use the conformal time r= j a _1 dt instead of the cosmic time t: 

D _ 1 D 

dt a dr 

D 2 _ 1 D 1 D a D ID 2 

d¥ ~ adr^adr' ~ ~a^d^ + ^ ' 



4 In general, the moduli space metric, h a b is neither canonically normalized nor field independent. How- 
ever, we show that in the class of models we are considering it is possible to make a field-dependent field 
redefinition such that the metric remains (approximately) flat throughout and after the inflationary period. 
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where ' denotes differentiation with respect to r and ' denotes differentiation with respect 
to t. Then, by making the change of variables, 5(f) 1 = ^itj, where the superscripts get 
changed to subscripts for later convenience, we have 

^- j ili u/d ui ui'a 
a a 2 a 2 a 2 

xll * o^ lh _L0 n/d2 Ul ' d UI " q h /i/o" 2 ' d \ fota\ 
dip 1 = 2—^- + 2 — 5 k- = —5 o—^ui + (2-^ ^-JM/ . (2.18) 

a or or a A or a A a a a A 

We also use the slow-roll approximation during inflation, 

e = ^«l (2.19) 
and integrate the conformal time by parts [16] [17], 

r ~ - ^ + 0(e 2 ) . (2.20) 
Thus, in terms of the conformal time (2.20), Eq. (2.16) becomes 

u'l + (k 2 - \)m = ^J2 MijUj ( 2 - 21 ) 
T J 

where the matrix Mjj is given by 

Mjj = eSjj + 2e u - VlJ - \tu{t- + + \ee u + 0(e 2 ) (2.22) 

cj) 1 cp J 

and the multi-field slow-roll parameters are defined as follows 

£/J -2l^ ( 2 - 23 ) 

= 1 d 2 V Vu 
mj Vdtfd^ V ' 

Note that the last term, |ee/j, is second order in slow-roll parameters and may be 
ignored. The forth term in (2.22) can also be treated as a second order term for the 
light fields. For the heavy fields, this term may be a first order term, ~ 0(l)ejj. For 
completeness sake, we will keep all the terms in our analysis throughout this paper. 

To solve (2.21), we have to decouple the equations by making a rotation U such that 

U~ 1 MU = diag{X I }, (2.24) 
where A/ are the eigenvalues of M. Then U is given by the similarity transformation 

U = ( g{ g~2 ■■■ g" n ) , (2.25) 



where gl are the eigenvectors of of M. 

Thus, by introducing the new fields, vi, 



ui = Uuvj , or vi = Ujjuj , (2.26) 
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we get that Eq. (2.21) is decoupled 



v >> + {k 2_f±A )vi = 0j (2 .27) 

where 

M/ 2 + 3A 7 . (2.28) 

Now we want a solution satisfying the Minkowski-like vacuum initial conditions [5] 
when kri 3> 1 {k 3> aH) 

uj =■ . (2.29) 

V2k 



The solution is _ 

vi = ^-^{-r^H^i-kr^m , (2.30) 

where ej are the normalized Gaussian variables 5 , satisfying [13] [18] [19] 

<ej(fc)) = (2.31) 
(ei(k)tj(k')) = 5u5 3 (k - k') . (2.32) 

2.3 The Asymptotic Solution 

We are mainly interested in the solution after Hubble exit when k < aH or kr — > 0. 
For small A/, 6 



fii = ^ + A/ + 0(A/ 2 ), for A/<1 



Eq. (2.30) becomes 



vi ~ ie^ Xl (l + C\ I )-^{-kT)- 1 ~ Xl e I {k), for kr -> , (2.33) 
V2A; 



where C = 2 — /o^2 — 7 (7 is the Euler-Mascheroni constant). 

For large and negative A/, 7 the order /i/ is complex (with a large imaginary part) and 
\fii\ » 1. We need to expand the Hankel function of large complex order. (For more 
details, see Appendix A). 

Recall (2.28), when A/ is large and negative, 



x l 9 - + 3Aj w ±ipj, p 7 > 0, » 1 , (2.34) 



5 The fluctuation can be treated as a random field which is a Gaussian process. The homogeneous 
universe can be divided into a set of sample space with different values of random fields mapped on it. 

6 These solutions are related to the light fields. Recalling (2.22), all the components in Mu related to 
the light fields are first order in slow-roll parameter. Thus the corresponding eigenvalues Xi are small, too. 

7 These solutions are related to the heavy fields. For the heavy fields, the relevant components in 
M/j(and thus Xi) are dominated by the diagonal elements of the matrix (Mu) which are the curvature of 
the potential, ~ ~^y L ^ 
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which is (almost) purely imaginary. Using (A-6), the solution (2.30) becomes 



where 8 



vi ~ -^(-r)5(l + pjy^e^-^P'e^^toe^k) (2.35) 
v2 



u _ e -p^^^-ip e i(^-p+plog\/i+p 2 ) _ ^yp e ~i(^-p+plog^i+p 2 ) (2.36) 



and z = —kr . 

Note that \cj\ ^ 1 in general. But for large p, 



\L0 



2 



1 + 2cosP(z)e- /m + e- 2pn « 1 (2.37) 



where (3{z) is some function of z = —kr according to (2.36). As we can see, the dependence 
on k for \u\ 2 (and hence the power spectrum), mainly given by the second term in (2.37), 
is exponentially suppressed. Hence, in the limit pi — > oo 

VI oc (_ r )|(l + /0 2)-| e l+(a-f ) PI PiJ^o ( _ r) l e 

and the solution for large p is suppressed by a factor of ~ 

This asymptotic solution can also be partially obtained from the following considera- 
tion. Consider the perturbation equation (2.27) when fij is large and negative and kr is 
small. It is approximate to the equation 

v'i + (-4>i = , 

z z 

where we change the variable r — > z = —kr. The solution of this equation is 

Vl(z) (X Z2 z V , 

which behaves similarly as (2.35). 

In summary, the perturbation solutions are 

r ie *f a j(1 + CAj ) i (_ fcr )-i-A Jej(A . )j )Aj | « ! ; 



Aj > 1 



where = ^/-(f + 3Aj) . 

3. The Curvature and Isocurvature Perturbations 
3.1 The Perturbations and Power Spectra 

It is convenient to decompose the scalar field perturbations into adiabatic perturbation 9 (parallel 
to the background trajectory) and entropy perturbation 10 (orthogoiaal to the background 

8 The dimension of (2.35) is |r| 5 . while the dimension of (2.33) is They are the same since kr is 

dimensionless. 

9 Also called curvature perturbations. 



10 



Also called non-adiabatic perturbations or isocurvature perturbations. 
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trajectory). We then define the adiabatic component 

5a = ^25(p I cos6 I (3.1) 

and the non-adiabatic component 

Ss 2 = J^j 2 -<&7 2 , (3.2) 
I 

with 



ol 91 91 (o o\ 

cos 6 = — ; = — (3.3) 

and 



>>'=-U,jv,. (3.4) 
a 

Note that in (3.4), the sum over the vj is for the light solutions only, since the perturbations 
of the heavy fields are strongly suppressed (see the discussion in Section 3.2 for more 
details). 

By definition, the two-point correlation functions (i.e., power spectra) are given by 

C xy 6 3 (k-k') = — 3 (x(k)y*(k')) (3.5) 
where x, y = 5a, 5s. For example, 

C aa (k)6 3 (k -k') = ^ (5a*5a) 

27T Z 



2ir 2 

1.1 

A; 3 

= ^ COS0 1 cos9 J U Ip U Jp (v*v p ) (3.6) 

/7T a Up 

where we used 

5(f)! = -Uijvj 
a 

For vj in the above expression, we have to omit the heavy solutions (where — Aj > 1) since 
they decay rapidly due to the expand of the universe (see discussions in the next section) . 
We next turn to the calculation of the power spectra using these correlation functions. 

3.1.1 Curvature Perturbations 

The comoving 11 curvature perturbation is defined by 

K = i) + —5a (3.7) 
a 



"Comoving means absent of peculiar motion. Comoving observers, such as large galaxies and galaxy 
clusters, measure zero momentum density at their own positions [4]. Their position, x, is time-independent 
in the unperturbed universe. Their physical coordinate is a(t)x. 
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In spatially flat gauge (tp = 0), the curvature perturbations become [5] [13] 

TZ=^5a. (3.8) 
a 

In multi-field inflation, in addition to the curvature perturbation, the isocurvature 
perturbations arise from the fluctuations orthogonal to the background trajectory 

S = —5s. (3.9) 
a 

The power spectrum of 1Z is defined as the expectation value of the Fourier components, 
which is just the ensemble average of the perturbations 

(n k K k ,)» = ^-PniWik - k% (3.10) 
Pn{k), = %C aa {k%. (3.11) 

Because of slow-roll approximation, the spectrum is usually calculated at Hubble crossing, 
denoted by *. In practice, Hubble crossing is often taken to be 50 or 60 e-foldings before 
the end of inflation[20, 21]. Due to the presence of isocurvature perturbation, the spectrum 
can change after Hubble crossing, which will be discussed in the following section. 
The power spectrum can be expanded around some ko [22] [23] 

Pn(k) = Pn(ko)(lr {ko) - 1+ " aln ^ , (3.12) 
ko 

where 

dlnP\{k) _ dn s 
n s {k) - 1 = — \ , a = — — . (3.13 
dink dink 

We have assumed that the momentum dependence of the running, a, can be neglected. In 
addition, a itself is of second order in slow-roll and should be small. We next turn to the 
power spectrum of the isocurvature fluctuation, P$, and the correlation power spectrum, 
Pus- 

3.1.2 The Isocurvature Perturbations 

The power spectrum of the isocurvature fluctuation, P$, and the correlation power spec- 
trum, Pus-, can be obtained in a similar way to the curvature perturbations. The non- 
adiabatic component has the general form 

5s = 5si = ^Sfa (3.1) 

For example, in a four-field model containing two heavy fields (0i, fo) and two light fields 

(03>4), 

(5 1 = (1,1, -cose 4 ,cos9 3 ) 



where 



cosO = — 
a 
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In a three field model, where there are two heavy fields (0i, <f>2) and one light fields ^3, 

(3 1 = (1,1,0) 

For the decoupled CclSG, clS will be discussed in section 3.2, we can totally ignore the 
heavy fields, and the coefficients reduce to the simpler forms 



P 1 



' '(0,0, -cos6 4 ,cos9 3 ) 
(0,0,0) 



for the four- and three-field models respectively. 
The perturbations are then given by 



5<f)i = -UjjVj, summed over the light vfs , (3-2) 
a 

from which the correlation functions follow 

(5s*5s) = (pSfitfSfa) = ^l3 J U Ip Uj p (v* p v p ) (3.3) 

(5a*5s) = (cosO 1 dtfp 3 'S<f>j) = ^cosO 1 [3 J U Ip U Jp (v*v p ) . (3.4) 

As before, the heavy v p s are ignored in the calculation. 

Thus, the power spectrum of the isocurvature fluctuation, P$, and the correlation 
power spectrum, Prs, at Hubble crossing are given by 

Ps(k)* = %C ss {k% (3.5) 

Cns(k)* = *^C as (k)U (3.6) 



where the two-point functions are given by 



C ss (k)5 3 (k - k') = —2 (5s*5s) = --^/3 t /3 j Ui p Uj p {v* p v p ) (3.7) 



k 3 k 3 

C as (k)5 3 (k -k') = ^ (5a*5s) = 1 ^cos0 I p J U Ip Uj p (v* p v p ) . (3.8) 

For future reference, it is convenient to define a dimensionless measure of the correlation 
angle between the power spectra [13] , 

cos A = Sjg S ; . (3.9) 
3.2 The Evolution of Perturbations After Hubble Exit 

For purely adiabatic perturbations, the curvature perturbation is a constant on super- 
horizon scales during the primordial era 12 [4] [5]. In this case, the observable perturbations 



12 The primordial era is defined as the period between Hubble exit and Hubble entry when the comoving 
scale, equals the Hubble scale, £ = A. 
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are calculated at horizon crossing. However, as Wands et al. have pointed out [13] [24] [25], 
the presence of entropy perturbations can change the curvature perturbation. In general, 
the time dependence of the curvature and isocurvature perturbation has the following 
form [25] [26] 

TZ = aHS (3.10) 



or in terms of the transfer functions 




S = I3HS (3.11) 



(3.12) 



The curvature perturbation on super-horizon scales is conserved if the perturbations are 
purely adiabatic or if the non-adiabatic perturbation is negligible. This general conclusion 
does not even depend on the slow-roll approximation or the form of the gravitational field 
equations (the specific theory of gravity) [24] . 

As we can see from the solutions of the perturbation equations (2.33) and (2.35), for 
each scale (1/fe), the spectrum of the perturbations with — Aj > 1 decay rapidly as the 
universe expands, ( The spectrum of the perturbations with |Aj| << 1, on 

the other hand, changes slowly, 

in the slow-roll parameters and the masses of the light fields over Hubble parameter. Thus 
we can ignore the contributions from the former and simplify the calculation. 
Recall (2.26) or 

5<f>i = -Uijvj 
a 

where Ujj is the transfer matrix determined by the mass matrix M of (2.22). If we 
assume that the heavy fields and the light fields are decoupled in such a way that the 
cross components Mjj (or Mjj), with / and J identified as light fields and heavy fields 
respectively, are subdominant compare to the non-cross components, then 

5<fti l « — C/jjjjWj,, Ii&cJi denote the light fields, (3.13) 



and 



-Uj h j h vj h , Ih&Jh denote the heavy fields. (3-14) 



This is true for most inflationary models encountered so far. For counterexamples, one has 
to use the full transfer matrix as in (2.26). Under the above assumption, the perturbations 
of the light fields, (|<50/ ; | 2 ) ~ (l^/J 2 }* [1 + 0(e) + O(^)], decay much slower than the 
perturbations of the heavy fields, {S<j>i h \ 2 ) ~ \- Therefore, in this case, one can neglect 
the contributions from the heavy fields when we calculate the curvature and isocurvature 
perturbations since the fluctuations in the heavy fields are strongly suppressed 13 . 



13 We can always do this unless the amplitude of the non-adiabatic fluctuation is greatly amplified at the 
end of inflation in the preheating stage[24][28][29]. 
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If there is a single light field (with all other fields being heavy) , the perturbations are 
purely adiabatic and the comoving curvature perturbation remains constant during infla- 
tion. If there is more than one light field, the cosmological inflation is driven by all the 
light fields. In addition to the adiabatic perturbation, they also produce entropy pertur- 
bation orthogonal to the background trajectory. In this case, the curvature perturbation 
is no longer a constant on super-horizon scales during inflation. The coupling between 
the entropy perturbation and the adiabatic perturbation, given by the (3.10) and (3.11), 
determines the evolution of the perturbations during and after inflation. 

In a typical two light field inflationary model, for example, with arbitrary potential 
and arbitrary background trajectory, it was shown [13] [25] that the scale-dependence of the 
observable spectra is determined by the slow-roll parameters at Hubble exit and the current 
observable cross-correlation. The amplitude of the power spectra are determined by the 
power spectra calculated at Hubble exit and the transfer functions which parameterize the 
detailed physics after Hubble exit until the end of reheating, given by (3.12) [25] 

P n = (1 + T% L s)Pn* + ZTnsCns* 
P S = Tl s P s * 

C-jzs = TssC-jis* + T-jisTssPs* (3.15) 
4. The Potential and Kahler Moduli Stabilization 

In what follows we focus on a particularly inflationary model derived from string theory 
consisting of multiple Kahler moduli, in the large volume limit (also known as the Large 
Volume Scenario) [7] [27]. We adopt the model originally proposed by Conlon and Quevedo 
in [9] and subsequently studied in [10] [11]. For more details, and in particular the conven- 
tions, see [27]. 

Supergravity in a four dimensional theory with M = 1 supersymmetry is completely 
specified by a Kahler potential fC and superpotential W. The Kahler potential is a real 
function of the complex scalar fields, while the superpotential is a holomorphic function 
depending only on the (ft 1 , and not the complex conjugate, (ft 1 . Focusing on the dynamics 
of the scalar fields relevant for inflation, the supergravity action is (we will work in the 
Einstein frame, and in units where Mp = 1) 



SV=i = / d^xyf^g 



(4.1) 



The scalar potential depends on the superpotential W, the Kahler potential K as well as 
the Kahler metric Qq, 

V((j>i, fr) = e K (g^DiWDjW - 3Ww) + V upMt (4.2) 
DiW = diW + WdiK (4.3) 
G fj = didjlC (4.4) 

The derivatives di and differentiate with respect to the 4n and 4q dependence, respec- 
tively. By expanding the complex fields in terms of their real and imaginary part, we can 
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relate the supergravity action above, (4.1), to the action, (2.1), discussed in section 2. The 
term Kiplift will include the effects of supersymmetry breaking arising from other sectors 
of the theory. 

We will demonstrate our methods in the context of Type IIB string theory compactified 
to four dimensions on a Calabi-Yau orientifold because the scalar potential in this case is 
well-understood and realistic four-dimensional models can be constructed [6, 7, 27, 30, 31, 
32]. After including the leading perturbative and non-perturbative corrections of string 
theory, the Kahler potential and superpotential are given by 

/C = -2 In [v + ^ j - ln(-i(r - f )) - In (^-i J Q A nj , 

I^GaAft + ^T^e-^ (4.5) 



W- 9 ' 



(4.6) 



Here g s is the string coupling, l s is the string length, 0, is the holomorphic three-form on 
the Calabi-Yau manifold M, G3 is the background field (flux) that is chosen to thread 
3-cycles in M and 

C(3) X (M ) 
? 2(2vr) 

where \ ls Euler number of M. The axion-dilaton field is r = Co + i e~^, and the integrals 
involving 0, are implicitly functions of the complex structure moduli. The fields Tj = 
Ti + ibi are the complexified Kahler moduli where t% is a 4-cycle volume (of the divisor 
Di € H^MjTj)) and bi is its axionic partner arising ultimately from the 4- form field. 
Here <2j = 2ir/Ni for some integer iVj, for each field, that is determined by the dynamical 
origin of the exponentials in the superpotential (iVj = 1 for brane instanton contributions, 
Ni > 1 for gaugino condensates). Finally, V is the dimensionless classical volume of the 
compactification manifold M (in Einstein frame, but measured in units of the string length). 
In terms of the Kahler class J = Yli^^i (by Poincare duality Di G H 2 (M, Z)), with the 
t l measuring the areas of 2-cycles, Cj, 

V = I j3 = !'••-.;' /;/,// . ( 4 - 7 ) 
Jm o 

where are the intersection numbers of the manifold. V should be understood as an 
implicit function of the complexified 4-cycle moduli via the relation 

n = d ti v = l -K^n k . (4.8) 

There are additional perturbative corrections to /C in (4.5), but we have kept the terms 
that give the leading contributions to the scalar potential in the large V limit of interest 
to us [33] . In particular, expanding fC to linear order in £ gives a consistent approximation 
in inverse powers of V. We have also assumed that all of the Kahler moduli Tj appear in 
the superpotential (see [31] for examples) and that we use a basis of 4-cycles such that 
the exponential terms in W take the form exp(-aiTi). As these exponentials arise from 
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an instanton expansion, in order to only keep the first term as we have done, the 4-cycle 
volumes must be sufficiently large to ensure that a{Fi S> 1. 

Finally, the form of the term VupMt in (4.2) depends on the kind of supersymmetry 
breaking effects that arise from other sectors of the theory. We take 

Uplift = ^2 ( 4 - 9 ) 

which will describe the energy of a space-filling antibrane [6] , fluxes of gauge fields living 
on D7-branes [34], or the F-term due to a non-supersymmetric solution for the complex 
structure/axion-dilaton moduli [35]. 

It was shown in [6] that a generic choice of background fields G3 causes all the complex 
structure moduli and the axion-dilaton to acquire string scale masses without breaking su- 
persymmetry. They are then decoupled from the low-energy theory and their contributions 
to K, and W are constants for our purposes 14 : 

K = -2M(v + £)-ln(£) +K o, 

W = 75 { W ° + ^ A ' e ~"" T ') ' (4 ' 10) 

where /Co (Wo) is the complex structure Kahler potential (superpotential), evaluated at the 
locations where the complex structure moduli have been fixed. It was shown in [7] that, 
when the Euler number, x < 0, for generic values of Wo (and hence of the background 
fluxes G3), the scalar potential for the Kahler moduli has a minimum where the volume V 
of the Calabi-Yau manifold M is very large - the associated energy scale is a few orders 
of magnitude lower than the GUT scale. Furthermore, in these Large Volume Scenarios 
there is a natural hierarchy - one of the Kahler moduli is much larger than the others and 
dominates the volume of the manifold. For our purposes they are also attractive because 
the scalar potential admits an expansion in inverse powers of the large volume V. This will 
allow us to carry out analytical calculations of inflation arising from Kahler moduli rolling 
towards the large volume minimum of the potential. 

Several previous works have considered inflation in the large volume setting, e.g., [9, 
36, 37, 38]. Here we include all Kahler moduli, and not just the light modes. Although 
we find that the heavy modes, corresponding to Kahler moduli that are stabilized before 
inflation takes, do not affect the dynamics during inflation in the models that we have 
studied, these modes do change after inflation has ended. 

Slow roll inflation can occur in a region of the field space where the potential is positive 
and very flat. We will look for this in the Large Volume Scenarios described above, where, 
at the minimum of the scalar potential, there is a hierarchy amongst the Kahler moduli 

Tl > T2,T 3 ,T 4 • • • (4.11) 



14 In the case of the F-term breaking due to the complex structure/axion-dilaton moduli [35], the con- 
tribution of the complex structure and axion-dilaton moduli to the scalar potential does depend on the 
volume (4.9). 
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which we will use to simplify the effective potential. 

For transparency of the equations, we will assume that the intersection numbers kijk 
are such that in the basis of 4-cycles, 73, the volume takes the diagonal form [27] 

V = a(n§ -£) W) = "«E A ^ § ( 412 ) 

i=2 t=l 

where Ai = —1, and \,i > 2 are usually positive. 

With the volume taking the above form we can explicitly compute the metric on the 
moduli space, Qg. = didjIC, which is needed both for the metric, h{j, and for the scalar 
potential, V, appearing in the four dimensional action (2.1). By expanding in inverse 
powers of V, keeping terms to 0(V~ 2 ), we obtain 

3aXi - 9a 2 AjAj 



g - = S 13 + - 3 \ 1 . (4.13) 

13 8(V+f)7i* 8(V+|)2 

With the axions minimized in the potential, the effective potential then becomes [9] 

^2 f= 4 . o t ,,/2 



3V Aia v 2 4V 3 V 2 ' 1 J 

i=2 i=2 

where we have assumed that /Co can be chosen such that the overall scale of the potential 
is simplified, i.e., overall factors of g s and 2tt are not present. Here we have expanded V 
to 0(V~ 3 ) to include the leading a'-corrections, , as well as the uplift term, The 

parameters in the potential can be chosen and tuned under certain constraints [36] [37] [39] . 

To determine the local minimum (vacuum) of the potential we need to solve the equa- 
tions 

dV 

(4 - 15 > 

While it is difficult to get the analytical results 15 , these equations can always be solved 
numerically. 

It is more convenient to work in the canonical frame, rather than the form taken by the 
supergravity metric in Eq. (4.13), since we have already solved the perturbation equations 
in the canonical frame 16 . Although it is difficult to find the exact transformations which 
can diagonalize the metric, we do find a canonical frame which is a good approximation as 
long as T\ 3> Tj, of which the field space transformations are 



(1 /3Ai(l + 3Ai) 1 / \ , « 

>! = 11 log(n) (4.16) 



4X W, i>2 (4.17) 



3ri 



15 Altough one can make approximations to solve the minimum equations analytically as in [9] and [10], 
it is desirable to solve them numerically. As we can show by numerical analysis, the analytical solutions 
after approximation will likely spoil the results. 

16 Note that comparing the kinetic energy terms in the actions (2.1) and (4.1), respectively, we find that 
hij = 2Qfj, with Qfj given in (4.13). 
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During, as well as after inflation, the metric, in terms of the above redefined fields, <t> 1 , 
remains canonically normalized, to leading order in inverse powers of the volume. Consid- 
ering that the original metric, Qq is a Kahler metric which neither is field independent nor 
diagonal, this result is somewhat surprising. 

5. Model Study 

In general, a multi-field inflationary model should contain both the heavy fields and the 
light fields 17 . To obtain inflation we choose the initial conditions such that the light fields 
are displaced away from the local minimum and the heavy fields are at the corresponding 
local minimum once the initial values of the light fields are chosen. We expect that the 
heavy fields will be frozen as the light fields approach the minimum. As we will see later in 
the numerical analysis, the light fields carry all the kinetic energy and are responsible for 
the creation of inflation. The number of e-foldings or the duration of inflation is determined 
by how far away the light fields are displaced from the minimum. The heavy fields will 
only begin to move and oscillate together with the light fields around the local minimum 
shortly after the end of inflation. 

In what follows, we will discuss two example models based on the discussion in the 
previous section. In both cases there are two heavy fields /moduli. The former has a 
single light field (inflaton) and the latter has two. By assigning appropriate values to 
the parameters in the effective potential, we solve the background equations of motion 
numerically. Next, we perform the field transformation (4.16), (4.17) to get the the kinetic 
energy in its canonical form. Then we use the perturbation solutions(light) to compute 
the curvature and isocurvature perturbations. Finally, we calculate the spectra and tilts 
at Hubble exit. Our models can be easily reduced or generalized. 

5.1 The Three-Field Model 

Let us construct an inflationary model with two heavy moduli, T\ and t 2 , and a light 
modulus, T3. This is essentially the Conlon-Quevedo model [9]. However, we do not 
assume that the initial values of the heavy moduli are the same as the final values after 
inflation. 

The parameters in the effective potential (4.14) are set to be 

1 2-7T 2tt 

a = — -,a 2 = ,a 3 = ,A 2 = 0.2, A3 = 0.002 

9V2 300' 100' 

Ai = -1, A 2 = 0.1, A 3 = 0.010, W = 500, f = 40, 7 = 9.75 x 10~ 6 
With these parameters, the local minimum is at 

T lmin = 62100.7, T 2min = 234.1, r 3min = 69.0202 

17 As has been shown in [11], the fields that are heavy (light) during inflation may become light (heavy) 
after inflation ends. So the heaviness (or the lightness) of a field is determined not only by the corresponding 
parameters, but to a large extent also by its position/value in the field space. 
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The initial conditions imposed are 18 
ri(0) = 76212.1, r 2 (0) = 246.99, r 3 (0) = 472, n(0) = r 2 (0) = 0, ,f 3 (0) = -7.13 x 10~ 19 
Obviously, all Tj(0)'s are quite different from Tj m j n 's. 





5 65.0 63.5 so ee.5 



64.5 65.0 65.5 66.0 66.5 



(a) 



(b) 



Figure 1: Inflation in the last few e- foldings, a) The slow-roll parameter e. b) The potential V. 





(b) 



Figure 2: Evolution of the heavy fields in the last few e-foldings. a) t\. b) r%. 

Apparently, the heavy moduli are frozen until the end of inflation, and Tfi na i — > T m i n . 
All the moduli began to oscillate after inflation ends. The light moduli oscillate much faster 
than the heavy moduli. The generated inflation lasts about Ntat = 65.2 e-foldings. We 
calculate the spectral index for the curvature perturbation at A* = 60 e-foldings (counted 
backward from the end of infaltion) before the end of inflation 



n K = n K, = 0.9563 (the running a = 0), e = 2.69 x 10 



-13 



Note that the spectral index is calculated at the Hubble exit (denoted by *). In the 
calculation (by (3.6)), we consider perturbations of both the heavy moduli and the light 
modulus. The same will be done in the next section. 

18 The attractor behavior of the evolution equations will ensure the same terminal velocity shortly after 
inflation begins regardless of how we choose the initial velocities [40] . 
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Figure 3: Evolution of the light field r 3 in the last few e-foldings. 
5.2 The Four-Field Models 

Symmetric Case First, we want to have two identical/symmetric light fields for which 
both the parameters and the initial conditions are the same. We expect to get a straight 
line background trajectory. The parameters in the effective potential (4.14) are set to be 

a = ' 02 = is ' as = m M = °- 2 ' As = °- 001 ' Aa = °- 001 

\ 1 = -l, A 2 = 0.1, A 3 = 0.005, A 4 = 0.005, W = 500, £ = 40, 7 = 9.75 x 10~ 6 
With these parameters, the local minimum is at 

r lmin = 62100.7, T2min = 234.1, r 3min = 69.0202, T 4min = 69.0202 

The initial conditions imposed are 

n(0) = 76212.1, r 2 (0) = 246.99, r 3 (0) = 472, r 3 (0) = 472, 

fi(0) = f 2 (0) = 0,f 3 (0) = f 4 (0) = -1.71 x 10" 19 




64.5 155.0 as j 66.0 ee.5 



(a) 




64.5 65.0 65.5 66.0 66.5 



(b) 



Figure 4: a) The slow-roll parameter e. b) The potential V. 
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(a) 



(b) 



Figure 5: Evolution of the heavy fields in the last few e-foldings. a) iv b) r 2 . 
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(a) 



(b) 



Figure 6: Evolution of the light fields in the last few e-foldings. a) T3. b) T4. 



Again, the heavy moduli remain frozen until after inflation ends. The generated infla- 
tion lasts about Ntot ~ 65.2. The spectral indices, computed as in the example above at 
N* = 60 e-foldings before the end of inflation, are given by 



nrc, = 0.9621 (a n = - 1.2 x 10~ 15 ), 
n 5 . = 0.9621 (a s = - 1.2 x 10~ 15 ), 
cosA* = -1.5 x 10~ 5 ,e=2.69 x 10~ 13 



(5.1) 
(5.2) 
(5.3) 



The spectral indices at Hubble crossing satisfies n-ji t = ns* because we have a symmetric 
model. The correlation angle cosA* is close to zero, and it is consistent with [13] in which 



cosA* ~ —2Cr) 



ICcosQ 1 /3 J V U \* = -1.98 x 10" 



where C ~ 0.7296. 



Nonsymmetric Case This time, we have two light fields of which the parameters are 
identical but the initial conditions are different. The parameters in the effective poten- 
tial (4.14) are set to be 



Q 



1 2tt 2tt t 

, a 2 = — , a 3 = — , A 2 = 0.2, A 3 = 0.001, A 4 = 0.001 



9^2' z 300 : 



100 
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\ 1 = -l, A 2 = 0.1, A 3 = 0.005, A 4 = 0.005, W = 500, £ = 40, 7 = 9.75 x 10" 

Tlmin = 62100.7, T2min = 234. 1 , T^min = 69.0202 

The initial conditions are 

Ti(0) = 76212.1, r 2 (0) = 246.99, r 3 (0) = 472, r 3 (0) = 492, 



fi(0) =r 2 (0) =0,f 3 (0) 



-1.72 x 10 



-19 



r 4 (0) 



-1.5 x 10 



-19 




64.5 65.0 



(a) 



(b) 



Figure 7: a) The slow-roll parameter e. b) The potential V. 
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65.0 65.5 66.0 ^-056.5 



(b) 



Figure 8: Evolution of the heavy fields in the last few e-foldings. a) T\. b) r%. 

The generated inflation lasts about N to t ~ 66. In this example, strictly speaking, 
inflation does not end initially when the innaton(r3) begins to oscillate. It ends when other 
fields also begin to oscillate. 

The spectral indices at N* = 60 e-foldings before the end of inflation are found to be 



= 0.9639(5^=4.69 x 10 -5 ), 
n St = 0.9879 (a.; =4. 85 x 10" 5 ), 
cosA* = -0.00501,e=1.44 x 10~ 13 . 

Again, the correlation angle cosA* is very small, consistent with [13] where 
cos A* ~ -2Cr] aSif = -2CCOS9 1 [i J V u U = 0.00507 



(5.4) 
(5.5) 
(5.6) 
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(a) 




64.5 65.0 65.5 66.0 66.5 



(b) 



Figure 9: Evolution of the light fields in the last few c-foldings. a) T3. b) T4. 



6. Conclusions 



In the typical inflation scenery, some or all of the light fields act as inflatons which are 
initially displaced from their local minimum. Along the path of the light fields the potential 
is very flat, and the light fields will slowly roll to their local potential minimum during 
inflation. Their perturbations are almost invariant and only decay slowly as the universe 
expands. 

The heavy fields, on the other hand, should be frozen during inflation, until shortly 
after the end of inflation when they start to move from the initial local minimum to the 
final local minimum. Their perturbations (as shown by the solution (2.35)) decay rapidly, 
approximately as a -3 , due to the expansion of the universe. That is why we ignore the 
heavy solutions (those with —A/ > 1) when we calculate the perturbations analytically. 

The resulting power spectrum of the curvature perturbation can be calculated us- 
ing (2.33), which is almost scale invariant due to the small A/. If there is only one light 
field, the contribution from the perturbations of the light fields is purely adiabatic. The 
perturbations are constants during inflation and reheating [5]. If there is more than one 
light field, they will also generate perturbations orthogonal to the background trajectory. 
In this case, the presence of the non-adiabatic/entropy perturbations can change the cur- 
vature perturbations, disqualifying it as a constant of time as discussed in section 3. 

The ratio of the curvature perturbation and isocurvature perturbation is 

Observational constraints from WMAP temperature anisotropy favor a small a with an 
upper limit 0.070(95% CL) [41]. The single (light) field inflation always satisfies this limit 
since there is no isocurvature perturbation. The multi-light-filed inflation can also satisfy 
this limit if the transfer functions can efficiently suppress the amplitude of the isocurvature 
perturbation after Hubble crossing. As discussed in section 3.2, the transfer functions 
depend on the details of the physics during as well as after the primordial era, and we hope 
to address this issue in the near future. 

We discussed an interesting inflation model derived from string compactification in the 
last two sections. Its potential has a nice shape which is ideal for generating inflation. The 
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result is very attractive. By adjusting the parameters and assigning appropriate initial con- 
ditions, we get the required number of e-foldings before the end of inflation. The numerical 
analysis shows that there is an oscillating period after the end of inflation which should 
be identified with the "preheating" era. It shows that the light fields oscillate much faster 
than the heavy fields. The preheating lasts a few e-foldings before the moduli roll to the 
minimum (vacuum) in the potential. In the single light field model, the spectral index cal- 
culated at iV* = 60 e-foldings before the end of inflation is about 0.956, which is consistent 
with the five-year WMAP data, n s = 0.960 ± 0.013 [21]. In the multi-light fields model, 
we get n s % ~ 0.962 — 0.964. However, in this case, it will require a better understanding 
of the transfer function, which is one of the main tasks left for future investigation, before 
we can compare the results with observations. 
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Appendix 

A. The Hankel Function of Large Complex Order 

First, we use the Frobenius expansion of the Bessel function around the origin 



■y.) = r 



1 1 /Z\2 1 

+ 



ro + i) r(^ + 2)V2; 2!r(^ + 3)V2 



To the lowest order, 



21 £ <-1)t f(T+7+l)t!' 2<<1 <^> 



■ / "-r(i + /i)" / -"-r(i-/i)/' z<< ( j 

So 

— ifiiv {z/2Y 0/2)"^ 

H Q)t z \ = e "t^Jjjjz) - J^{z) _ e rg+g r(i- M ) (A _ 3) 
^ — isimifi —isimrfi 

Stirling's formula for large v approximation of the Gamma function gives 



2vr r . 1 1 



T{u) = v v e-^ — [1 + — + _ + 0(1,-3)], \ v \ >l, \arg(u)\<7r-e (A-4) 
where e is any small positive number. 
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Let v = 1 ± ip where p > 1 , to lowest order, 




= (1 ± ip)(l ± ip^Pe^eVy/tor 



1 



(l + p 2 )i 



\/2tt(1 + p 2 )Je 



l-ap e ±i(f-p+pio 9 ^/T+^) 



(A-5) 



where 1 ± ip = -y/l + p 2 e : 



So for p > 1, the Hankel function is 




f(l+p 2 )-3 e 1+a V 



H = ip 



H = -ip 



(A-6) 



where 




(A-7) 



References 

[1] A. H. Guth, "The Inflationary Universe: A Possible Solution To The Horizon And Flatness 
Problems," Phys. Rev. D 23, 347 (1981). 

[2] A. D. Linde, "A New Inflationary Universe Scenario: A Possible Solution Of The Horizon, 
Flatness, Homogeneity, Isotropy And Primordial Monopole Problems," Phys. Lett. B 108, 
389 (1982). 

[3] A. J. Albrccht and P. J. Stcinhardt, "Cosmology For Grand Unified Theories With 
Radiatively Induced Symmetry Breaking," Phys. Rev. Lett. 48, 1220 (1982). 

[4] D. H. Lyth and A. Riotto, "Particle physics models of inflation and the cosmological density 
perturbation," Phys. Rcpt. 314, 1 (1999) [arXiv:hep-ph/9807278]. 

[5] D. Baumann, "TASI Lectures on Inflation," arXiv:0907.5424 [hep-th]. 

[6] S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, "De Sitter vacua in string theory," Phys. 
Rev. D 68, 046005 (2003) [arXiv:hep-th/0301240]. 

[7] V. Balasubramanian, P. Bcrglund, J. P. Conlon and F. Quevedo, "Systematics of Moduli 
Stabilisation in Calabi-Yau Flux Compactifications," JHEP 0503, 007 (2005) 
[arXiv:hep-th/0502058]. 

[8] V. Balasubramanian and P. Bcrglund, "Stringy corrections to Kahler potentials, SUSY 
breaking, and the cosmological constant problem," JHEP 0411, 085 (2004) 
[arXiv:hcp-th/0408054]. 

[9] J. P. Conlon and F. Quevedo, "Kachler moduli inflation," JHEP 0601, 146 (2006) 
[arXiv:hep-th/0509012]. 

[10] H. X. Yang and H. L. Ma, "Two-field Kahler moduli inflation on large volume moduli 
stabilization," JCAP 0808, 024 (2008) [arXiv:0804.3653 [hep-th]]. 

[11] J. J. Blanco-Pillado, D. Buck, E. J. Copeland, M. Gomez-Reino and N. J. Nuncs, 
arXiv:0906.3711 [hep-th]. 



- 23 - 



[12] C. P. Burgess, J. M. Clinc, F. Lemieux and R. Holman, "Are inflationary predictions 
sensitive to very high energy physics?," JHEP 0302, 048 (2003) [arXiv:hep-th/0210233]. 

[13] C. T. Byrnes and D. Wands, "Curvature and isocurvature perturbations from two-field 

inflation in a slow-roll expansion," Phys. Rev. D 74, 043529 (2006) [arXiv:astroph/0605679] . 

[14] V. F. Mukhanov, H. A. Feldman and R. H. Brandcnberger, "Theory of cosmological 

perturbations. Part 1. Classical perturbations. Part 2. Quantum theory of perturbations. 
Part 3. Extensions," Phys. Rcpt. 215, 203 (1992). 

[15] T. T. Nakamura and E. D. Stewart, "The spectrum of cosmological perturbations produced 
by a multi-component inflaton to second order in the slow-roll approximation," Phys. Lett. B 
381, 413 (1996) [arXiv:astro-ph/9604103]. 

[16] D. J. Schwarz, C. A. Terrero-Escalantc and A. A. Garcia, "Higher order corrections to 
primordial spectra from cosmological inflation," Phys. Lett. B 517, 243 (2001) 
[arXiv:astro-ph/0106020] . 

[17] S. Hirai and T. Takami, "Scale dependence of the power spectrum of the curvature 

perturbation determined using a numerical method in slow-roll inflation," arXiv:0908.0585 
[astro-ph.CO]. 

[18] V. F. Mukhanov, Physical Foundations of Cosmology, Cambridge University Press, 2005. 

[19] Z. Lalak, D. Langlois, S. Pokorski and K. Turzynski, "Curvature and isocurvature 

perturbations in two-field inflation," JCAP 0707, 014 (2007) [arXiv:0704.0212 [hep-th]]. 

[20] E. F. Bunn, A. R. Liddle and M. J. . White, "Four-year COBE normalization of inflationary 
cosmologies," Phys. Rev. D 54, 5917 (1996) [arXiv:astro-ph/9607038]. 

[21] E. Komatsu et al. [WMAP Collaboration], "Five- Year Wilkinson Microwave Anisotropy 
Probe (WMAP)Observations:Cosmological Interpretation," Astrophys. J. Suppl. 180, 330 
(2009) [arXiv:0803.0547 [astro-ph]]. 

[22] S. Hirai and T. Takami, arXiv:0812.1079 [astro-ph]. 

[23] H. V. Peiris et al. [WMAP Collaboration], "First year Wilkinson Microwave Anisotropy 
Probe (WMAP) observations: Implications for inflation," Astrophys. J. Suppl. 148, 213 
(2003) [arXiv:astro-ph/0302225]. 

[24] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, "A new approach to the evolution of 
cosmological perturbations on large scales," Phys. Rev. D 62, 043527 (2000) 
[arXiv:astro-ph/0003278] . 

[25] D. Wands, N. Bartolo, S. Matarrese and A. Riotto, "An observational test of two-field 
inflation," Phys. Rev. D 66, 043520 (2002) [arXiv:astro-ph/0205253]. 

[26] N. Bartolo, S. Matarrese and A. Riotto, "Adiabatic and isocurvature perturbations from 
inflation: Power spectra and consistency relations," Phys. Rev. D 64, 123504 (2001) 
[arXiv:astro-ph/0107502] . 

[27] J. P. Conlon, F. Quevedo and K. Suruliz, "Large- volume flux compactifications: Moduli 
spectrum and D3/D7 soft supersymmetry breaking," JHEP 0508, 007 (2005) 
[arXiv:hcp-th/0505076]. 

[28] L. Kofman, A. D. Lindc and A. A. Starobinsky, "Towards the theory of reheating after 
inflation," Phys. Rev. D 56, 3258 (1997) [arXiv:hcp-ph/9704452]. 



- 24 - 



[29] C. T. Byrnes, "Constraints on generating the primordial curvature perturbation and 
non-Gaussianity from instant preheating," JCAP 0901, Oil (2009) [arXiv:0810.3913 
[astro-ph]]. 

[30] S. B. Giddings, S. Kachru and J. Polchinski, "Hierarchies from fluxes in string 
compactifications," Phys. Rev. D 66, 106006 (2002) [arXiv:hep-th/0 105097]. 

[31] F. Dencf, M. R. Douglas, B. Florea, A. Grassi and S. Kachru, "Fixing all moduli in a simple 
F-thcory compactification," Adv. Thcor. Math. Phys. 9, 861 (2005) [arXiv:hcp-th/0503124]. 

[32] J. P. Conlon, C. H. Kom, K. Suruliz, B. C. Allanach and F. Qucvedo, "Sparticle Spectra and 
LHC Signatures for Large Volume String Compactifications," arXiv:0704.3403 [hcp-ph]. 

[33] M. Berg, M. Haack and B. Kors, "String loop corrections to Kachlcr potentials in 
oricntifolds," JHEP 0511, 030 (2005) [arXiv:hep-th/0508043]. 

[34] C. P. Burgess, R. Kallosh and F. Quevedo, "de Sitter string vacua from supersymmetric 
D-terms," JHEP 0310, 056 (2003) [arXiv:hep-th/0309187]. 

[35] A. Saltman and E. Silverstein, "The scaling of the no-scale potential and de Sitter model 
building," JHEP 0411, 066 (2004) [arXiv:hep-th/0402135]. 

[36] J. R. Bond, L. Kofman, S. Prokushkin and P. M. Vaudrevange, "Roulette inflation with 
Kachlcr moduli and their axions," Phys. Rev. D 75, 123511 (2007) [arXiv:hcp-th/0612197]. 

[37] J. Simon, R. Jimenez, L. Verde, P. Berglund and V. Balasubramanian, "Using cosmology to 
constrain the topology of hidden dimensions," arXiv:astro-ph/0605371. 
V. Balasubramanian, P. Berglund, R. Jimenez, J. Simon and L. Verde, "Topology from 
Cosmology," JHEP 0806, 025 (2008) [arXiv:0712.1815 [hep-th]]. 

[38] R. Holman and J. A. Hutasoit, "Systcmatics of moduli stabilization, inflationary dynamics 
and power spectrum," JHEP 0608, 053 (2006) [arXiv:hep-th/0606089]. 

R. Holman and J. A. Hutasoit, "Axionic inflation from large volume flux compactifications," 
arXiv:hcp-th/0603246. 

[39] A. C. Vincent and J. M. Cline, "Curvature Spectra and Nongaussianities in the Roulette 
Inflation Model," JHEP 0810, 093 (2008) [arXiv:0809.2982 [astro-ph]]. 

[40] C. Ringeval, "The exact numerical treatment of inflationary models," Lect. Notes Phys. 738, 
243 (2008) [arXiv:astro-ph/0703486]. 

[41] C. Hikage, K. Koyama, T. Matsubara, T. Takahashi and M. Yamaguchi, "Limits on 

Isocurvature Perturbations from Non-Gaussianity in WMAP Temperature Anisotropy," Mon. 
Not. Roy. Astron. Soc. 398, 2188 (2009) [arXiv:0812.3500 [astro-ph]]. 



- 25 - 



